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Abstract 

The supersymmetry (SUSY) self-consistent approximation for the model of non-equilib- 
rium thermodynamic system with quenched disorder is derived from the dynamical action 
calculated by means of generalized second Legendre transformation technique. The equa- 
tions for adiabatic and isothermal susceptibilities, memory and field induced parameters are 
obtained on the basis of asymptotic analysis of dynamical Dyson equations. It is shown that 
the marginal stability condition that defines the critical point is governed by fluctuations 
violating fluctuation-dissipation theorem (FDT). The temperature of ergodicity breaking 
transition is calculated as a function of quenched disorder intensities. Transformation of 
supernelds related to the mapping between an instanton process and the corresponding 
causal solution is discussed. 
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1 Introduction 

In recent years significant study has been given to the microscopic theory of non- 
equilibrium thermodynamic systems with quenched disorder that reveal non-ergodic 
behavior and exhibit memory effects. Spin glasses [1,2] and random heteropoly- 
mers [3], that received the most of attention, provide the well known examples of 
such systems. Procedure of the averaging over disorder is at the heart of theoretical 
approaches developed for the description of the systems. 
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Thus the bulk of static theories are based on the replica method firstly introduced 
in [4] or employ the methods of random matrix theory (RMT) (see [5] for review). 
In addition, the powerful supersymmetry approach by Efetov [6], where generating 
functional describing statistics of the density of states is mapped on the nonlinear 
supermatrix cr-model, has had RMT as a key ingredient of its development. 

The starting point of dynamical theories is stochastic dynamics defined in terms 
of the corresponding stochastic equations. A wide range of problems can be for- 
mulated in such a way: kinetics of ordering [7], non-equilibrium dynamics of 
spin-glass systems [2,8,9] and so on. When the stochastic dynamics is governed 
by Langevin equations, the generating functional of the stochastic problem can be 
written as a functional integral [10,1 1] (MRS formalism) thus allowing for the av- 
eraging over disorder at the very beginning of calculations. It has led to the field 
theoretic formulation of the problem, so one could take advantage of using the 
machinery of the field theory. The theory is known can represented in the super- 
symmetric form that reveals hidden supersymmetry (SUSY) of the stochastic prob- 
lem [12,22]. 

In this paper SUSY formalism in the superfield representation is employed to study 
ergodicity breaking transition in the model thermodynamic system with quenched 
disorder. It implies that generating functional of Langevin dynamical system is rep- 
resented as functional integral over superfields with Euclidean action by means of 
introducing Grassmann anticommuting variables. These variables and their prod- 
ucts serve as a basis for superfields which components include Grassmann fields in 
addition to the usual real-valued fields. The correlation functions of the superfields 
(supercorrelators) then encode physically relevant information on observables in 
components of the correlators that are the autocorrelator of the order parameter 
field C and the response functions G± (advanced and retarded Green functions). 
The analysis is based on a system of dynamical Dyson equations supplemented with 
the equations for averaged order parameter provided the distribution of quenched 
random variables is not symmetric. It was pointed out in [14] that, in the mode- 
coupling approximation, the Dyson equations can be regarded as Euler-Lagrange 
equations for the functional of supercorrelators known as dynamical action and this 
functional bears striking similarity with the replica expression for the free energy. 

In our study the dynamical action is shown to be a second Legendre transforma- 
tion [15] of the free energy. So, it depends on both averaged order parameter and 
two-time correlation functions and can be calculated on the basis of suitably de- 
fined diagrammatics. On the technical side, since algebraic structure of the super- 
field representation is directly related to the underlying symmetry, we have the 
reduction of the number of diagrams to be taken into consideration in the perturba- 
tion theory. From the other hand, SUSY formalism allows ergodicity breaking tran- 
sition be interpreted as a dynamical symmetry breaking transition. Indeed, in the 
SUSY language, the well-known "causality" condition and fluctuation-dissipation 
theorem (FDT) immediately follow from the corresponding Ward identities [9,16]. 
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Below the critical temperature FDT and time translational invariance are dynam- 
ically violated that result in the appearance of anomalous solutions to the Dyson 
equations. In particular, the latter include the case, where the system being in a 
non-ergodic state is characterized by a very slow relaxation. This phenomenon is 
known as aging and, as it has been demonstrated for a number of solvable mod- 
els [17-19], the effect is a consequence of trapping in metastable attractors. 

In this article, leaving aside detailed study of aging, we shall perform asymptotic 
and stability analysis of dynamical Dyson equations to characterize the transition 
in terms of asymptotic quantities such as adiabatic and isothermal susceptibilities, 
memory and field induced parameters. This approach is applied to the simple model 
of thermodynamic system with quenched external field and two body interaction. 
Note that the second Legendre transformation technique plays the unifying role in 
this theory, so that the SUSY based theory can be directly employed for the study 
of ergodicity breaking transitions in different disordered systems such as random 
heteropolymers [3] and filled nematic liquid crystals [20]. 

Layout of the paper is as follows. 

In Sec. 2 the formalism of SUSY approach is briefly outlined. Second Legendre 
transformation for disordered systems is introduced in Sec. 3. In Sec. 4 the model 
of non-equilibrium thermodynamic system with quenched disorder is studied on 
the basis of asymptotic analysis of the Dyson equations. In high temperature region 
it is found that the relevant parameters are the static susceptibility x an d the field 
induced parameter q h . The latter is due to the presence of random field that affects 
asymptotics of the autocorrelator C(t). Equations for \ an d qh combined with the 
marginal stability condition define the temperature of ergodicity breaking T c . De- 
pendencies of T c on the quenched disorder intensities are calculated. It is shown 
that the low temperature region T < T c can be described in terms of adiabatic 
Xa an d isothermal \ susceptibilities, the dynamical Edwards-Anderson memory 
parameter q and the field induced parameter q h , so that the role of an order (non- 
ergodicity) parameter plays the difference Aq = q — q h . Discussion of numerical 
results and concluding remarks are given in Sec. 5. Appendix A details the remark 
that the superfield representation induced by the shift in time, t — > t — 66, leads to 
the mapping between an instanton process and normal downhill motion. 



2 General SUSY formalism 



In this section we sketch the general formalism of SUSY based theory of a non- 
equilibrium thermodynamic system. It serves as an introductionary part of the paper 
and gives some details on the results used in subsequent sections. For definiteness, 
in what follows we use the lattice designations, so that the field rji(t) defined on 
sites of the lattice (the sites are labelled with the index i) gives configuration of 
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order parameter at the instant of time t. 

Relaxational dynamics of the order-parameter field in the presence of thermal noise 
is governed by the Langevin equation: 

5V 

rn = -T- + Ci{t), CD 
or]i 

where the relaxation constant is absorbed by suitable rescaling of time and (. The 
thermodynamic potential V is assumed to be a t-local functional 

V{rj} = J dtV( V ), (2) 

and d(t) are Gaussian stochastic functions subjected to the white noise conditions: 

(C i (0C-(0)> C = 2TV(0» (0)c = 0, (3) 
where T is the temperature. 

2.1 Generating functional in superfield representation 

The starting point of the MRS formalism [10,11] is the generating functional for 
correlation functions of the stochastic problem written in the form: 

^w = (/n^ det (^) w^)) ex p(/ d ^)) » (4) 

where S(-) is the delta function, L(rf) = —d t r] — 5V/5r/ + ( and (-)^ denotes aver- 
aging over noise realizations. (Hereafter the functional notations will be adopted 
assuming summation over repeated indexes and integration over repeated non- 
discrete arguments.) 

Exponentiating the delta function through 'Lagrange multipliers' <fi{t) and the Ja- 
cobian functional determinant through Grassmann fields (ghosts) tpi(t), ^(t) and 
averaging away the noise we obtain functional integral representation for the gen- 
erating functional (see [21] for recent discussion of the corresponding steps) 

Z M = jl[Dr ]l Dip l Dtlj i Dtlj l exp{-S + J dtu lVl } (5) 
where the action reads 

S = J dtL, (6) 
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nn 2 f SV) - 

L = -Tw + <pi yii + — | - ipi 



d S 2 V 
+ 



dt 5i]i5r]j 



^3 



(7) 



Introducing the superfields 



4>i{z) = ^ = r}i + + + 99<p h z = {t, e, 9} = {t, 0}, 



(8) 



where 9 and 9 are anticommuting Grassmann variables, and substituting {9, 9} — > 
{T- 1 ' 2 9,T- 1 / 2 9} (or, alternatively, -> {T 1 / 2 ^, T 1 / 2 ^}, ^ -> ) we 

derive the action as a functional of superfields 



s = ^S dzL = fI dtd2 ° L > 



(9) 



L = D<j> i D<l> i + V{<j>), 



(10) 



where d 2 = dtfdtfand 



- d d d 

= W = d9~~dt 



(11) 



It is not difficult to show that the operators D and D enjoy the following properties: 



D 2 = D 2 = 0, {D,D} = -1 [D,Df = ^ 



(12) 



where the curly brackets stand for anticommutator, and the kinetic term in the action 
dzD(f)iD(j)i can be written in the form 



/ 



Note that the form of superfield representation (8) is fixed up to the change of the 
basis in the space of superfields, so that the corresponding transformation would 
give another representation of the SUSY group. Interestingly, as it is shown in 
Appendix A, alternative superfield representation, generated by the transformation: 
<f)(t, 0) — > (j)(t — 99, 0), can be used to construct the mapping between an instanton 
process and the corresponding causal solution inverted in time. 
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2.2 Symmetries and Ward identities for 2-point functions 



The action by Eqs. (9,10) is invariant under the action of the group of supersym- 
metry [9,14] with generators given by 

d d - d d 

D= W D= m +e 8i- W (13) 

As a consequence, the correlation function G(zi, z 2 ) = {(j>(zi)(f>(z2)) (for brevity, 
the indexes of superfields are suppressed) meets a set of Ward identities provided 
that the supersymmetry is not broken [9]. Clearly, the invariance under translations 
in time implies that G(z 1: z 2 ) depends only on t = ti — t 2 . Another two identities 
are of special interest to us: 

(D' 1 + D' 2 )G(z 1 ,z 2 ) = 0, (14) 



(D' 1 + D' 2 )G(z 1 ,z 2 ) = 0, (15) 

where 

f)'-— n>-— + f) — 

Eq. (14), known as "causality condition", implies that the correlator is of the fol- 
lowing form 

G( Zl ,z 2 ) = C(h,t 2 ) + (h - 9 2 ) {G + {hM)0i ~ G4h,t 2 )e 2 ) , (16) 

where 



C(t 1 ,t 2 ) = {r j (t 1 )r j (h)), (17a) 
G + (tx, t 2 ) = (tpihHh)) = m^ih)) , (17b) 
G-(ti,h) = (r ] (t l Mt 2 )} = -(^( tl )^(t 2 )} . (17c) 

Thus the identity (14) allows the correlators of Grassmann fields ip and ip be ex- 
pressed in terms of response functions. 

Algebraic structure of the correlator can be conveniently emphasized by introduc- 
ing a set of operators {T, A + , A_} 



A + (O,0') = (6-6')6, 

A_(0,0')=A + (0',0) = -( 
T(O,0') = l, 



- o'W, 



(18a) 
(18b) 
(18c) 
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equipped with the operator product 

A 1 -A 2 (e,e') = I d 2 e"A 1 (e,e")A 2 (e",e'), (19) 

so that 

G = CT + G+A+ + G_A_ . (20) 



It is not difficult to verify that the above operators form the basis of algebra with 
respect to the operator product 



A± • A± = A±, T A + = A T = T, (21a) 
A± • A T = A + • T = T ■ A = T ■ T = 0. (21b) 

In particular, Eqs. (21) ease finding of inversion formulae for superoperators. For 
example, we can derive the expression for G -1 : 

G 1 = G+ 1 A + + GZ 1 A_ - GZ 1 ■ C ■ G^T , (22) 

so that 

G G 1 = 5(6 - 6') = (0 - 9'){9 - 6') = A + + A_ . (23) 
Note that expansion of the bare correlation function G^ over the above basis is 

G (o) = p(2) + m | 5(6> _ y.. = 

= {—2 T + (m — d t ) A + + (m + d t ) A_}%. (24) 



The second identity Eq. (15) provides the relation between the autocorrelator of 
the order parameter C(t,t') and the response functions G-(t,t') = G(t,t') and 
G + (t, t') = G(t' : t) (retarded and advanced Green functions) known as fluctuation- 
dissipation theorem (FDT): 



SLc(t, t') = -6{t - t')G(t, t') + 6{t' - t)G(t', t) (25a) 

C6 

r) 

— C(t, t') = 6(t - t')G(t, t') - 6(t' - t)G(t', t) (25b) 
where 9(t) is the step function. 
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It is convenient to reformulate FDT (25) in terms of time-dependent susceptibility 

x(M') 

X (t,t')= fdr'G{t,r') (26) 
Jt' 

that gives the response to an applied field held constant from time t' up to t. (Here- 
after it is assumed that t > t' .) Since FDT implies the time translational invari- 
ance, G, C and \ depend only on the time separation r = t — t'. Then integrating 
Eq. (25b) and taking the limit r — > oo yields the required form of FDT: 

X = / dr'G(r') = q -q h , (27) 
Jo 

where x is the static susceptibility, q = C(t, t) = C(0) and q h = C(oo). 

Analogously, Ward identities for proper vertices lead to the same relations for the 
mass operator (self-energy), 

X( Zl ,z 2 ) = E(t 1 ,t 2 )T + E + (t 1 ,t 2 )A + + E_(t 1 ,t 2 )A_, (28) 

that enter the Dyson equation. 



3 Dynamical action for system with quenched disorder: second Legendre 
transformation technique 

In this section we derive the dynamical action as a second Legendre transformation 
of the free energy functional. Since it is straightforward to generalize the subse- 
quent considerations, for the sake of simplicity, the technique will be employed to 
study the system with the thermodynamic potential V(<f>) of the following form 

V(<t>) = £{Wi(*)) + HiMz)} - £ Wi^(z)^-(z), (29a) 

i ij 

U{<P) = + tW0) = |0 2 + , (29b) 

where m = fiT. Under suitable assumptions, it can be regarded as a discrete version 
of the well-known Landau-Ginzburg free energy functional for coarse-grained 
scalar order parameter field [7]. We consider the case where the couplings Wij 
and the field correspond to quenched degrees of freedom and are independent 
Gaussian variables: 

W^mn 0> Wmn^ m'n' ^mrn'^nn'U^mni H m 0, H m H m i 5 mrn 'h . 

Averaging away the quenched variables gives 
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exp{-S} = ex.-p{-T~ 1 J dzdz'L eff } (30a) 
L eff = E {\Mz)K(z, z'^z') + 5(z - z')U anh {Uz))} ~ 

-(2T)- 1 ]T u>ij<f>i(z)<f>j (z)<f>MMz') ( 30b ) 



where (see Eq. (24)) 



K(z,z') = 



2 



~T + (2T + (m - $)A+ + (m + d t )A_) <S(* - t') 



V (3D 



Note that, by contrast with the real two-particle interaction, the term generated by 
the averaging is non-local in time. 

At this stage we can use functional methods of the field theory to derive equations 
of motion for the one- and two-point correlators, (<j>i(z)) and ((pi(z)(pj(z')). To this 
end let us write the effective action (30) in slightly generalized form: 



1 ( 

s eff = Y^S ((l>i\AW) + -J dz 1 dz 2 ^w ij (t)i(z 1 )(f)i(z2)(f)j(z 1 )(l) j (z2) , (32) 



where 



5 (^|A«) = ^ - fdz! . . . dz n A^( Zl , . . . z n )Uzi) ■ ■ ■ Mzn) (33) 

71=1 U - J 



and A® = {A®, A%\ A®, A®} stands for a set of "superpotentials" that define 
the one-particle action So and are assumed to be symmetric functions of argu- 
ments. In our case A? = Af = 0, Af = -K(z u z 2 )/T and Af = -\/T5{z 1 - 
z 2 )S(z 2 — z 3 )5(z3 — z 4 ). For convenience, we will keep the generalized denotions 
and the above A$ will be inserted into the final equations. Note that A 1 is com- 
monly referred to as an external field and the inverse of operator with the kernel 
—A 2 (zi : z 2) is proportional to the bare correlation function. 



The effective action S e ff defines the partition function Z(A) and the free energy 
F(A) as generating functionals of n-point correlators without vacuum loops and 
connected n-point correlators, respectively: 
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Z(A) =eMF(A)} = N- 1 J J] ex P {S e// (0|A)} , 



8 n Z(A) 



5Af{ Zl )...8At\z, 
5 n F(A) 



(0/ 



<L4S l) (^)...cL4^(z 
<5F(A) 



(0/ 



*4°(*i, ^2) 



(&(zi) • • • = n\ o${z± ...z n ), 

(&(zi) ■ ■ ■ 4>i{z n ))c = (3 i n ) { Z l ■■■Zn) , 

a { 2 ) (z 1 ,z 2 ) ■ 



(34a) 
(34b) 
(34c) 
(34d) 



After introducing auxiliary field Xi(zi,z 2 ) and performing the Hubbard-Stratonovich 
transformation 



exp 



- J dz 1 dz 2 J2 W ij ( Pi(Zl)Mz2)(pj(z 1 )(f) j (z2) 
ij 



/ YiDXiCxp -\J dzidz 2 Y^ w i] lx i{zi, z 2 )Xj(zi, z 2 ) + 

i 1 ij 



+ X i ( Z l , Z 2 ) <Pi (Zl ) (pi (Z 2 ) 

i 

the partition function assumes the following form 



(35) 



Z ( A ) = J Il DX * ex P -\J dz 1 dz 2 J2^ 1 Mzi, Z2)Xj(z u z 2 ) + 

i L ij 



+ Y,Fo(4\A? + 2X l ,A?,A?) 



(0 4 (0> 



(36) 



where F (A^) is the generating functional of connected correlators for the one- 
particle action S : 



exp{F (A«)} = / Z^ex P {S (^|A«)}. 



(37) 



After changing variables A^ 2 = A^ 2 + 2X; in Eq. (36) it is ready to derive the 
equations of the saddle-point approximation: 



SFp 

5Aip 



(») 



(38) 



It is known that the mean field approximation is given by the leading order in steep- 
est descent calculations of such kind [16]. The only difference from the standard 
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mean field theory is that the interaction term in Eq. (36) contains A$ instead of 
A±\ So it is natural to define second Legendre transformation for F as a general- 
ized Legendre transformation with respect to the first two "superpotentials" A 1 and 
A 2 [15] : 



r(«i, «2|^, A 4 )=F(A) - Af^ - 4*) __ 



= F(A)-i4iai-A 2 a 2 , (39) 

where 

af { Zl ) = (3? { Zl ) = 2a%\z 1 ,z 1 )=$\z 1 ,z 2 )+(3?(z 1 )P?(z 2 ) 
and universal functional notations used in (39) imply that 

Aiai =YlJ dzxA^ (zi)a[ l) (zx) , A 2 a 2 = YlJ dz 1 dz 2 A i 2 ) (zi, z 2 )a%\z 1 , z 2 ) . 

i i 

Similar to the standard Legendre transformation, given the functional V with A 1 
and A 2 , the correlators a x and a 2 are defined as solutions to the equations 



' !l ~4\ Zl ) (40a) 



Sa?( Zl ) 



' )] -A?( Zl ,z 2 ). (40b) 



5a 2 \zi, z 2 ) 

Eqs. (40) are Euler-Lagrange equations for the functional 

$(ai, a 2 \A) = r(ai, a 2 |A 3 , A4) + aiAi + a 2 A 2 , (41) 

so that substituting of extremals from Eqs. (40) into $ provides the corresponding 
values of the free energy F. It follows that $ has the meaning of dynamical action. 

In the mean field approximation, from Eq. (38) the expression for T is 

T = J2^o(a?,a ( 2 l) \A 3 ,A 4 ) + 2j2w l3 J Az^af (z u z 2 )a^ (z u z 2 ) , (42) 

i ij 

where T is the Legendre transform of the one-particle free energy. As a functional 
of Pi and f3 2 , T is given by [15] 

r = iTrln/? 2 + f . (43) 

where f can be calculated as a sum of two-particle irreducible diagrams (it cannot 
be made disconnected by cutting off two lines) with "dressed" internal lines (f3 2 is 
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the propagator). Up to the second order of the perturbation theory with A 3 = 0, we 
have 

f o = ^A, (# + Q(3 2 (3l + 3Pl) + ^A 4 p 4 2 A 4 + l^O/^A^). (44) 
The diagrammatic representation of the terms on the right hand side of Eq. (44) is 



1 

24 




+ 
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+ 



+ 



48 
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In the diagrams zigzag lines correspond to /^i-lines directly joined to the vertices. 
The equations of motion in new variables are 



" ! 2 ^ 6T : fr{z 2 ) = -A 1 (z 1 ) 



5(5i(zi) 5(3 2 (zi,z 2 ) 

5p 2 {zi,z 2 ) 
Note that, owing to the identity 

5 



Trln/5 2 = (3 2 - 1 {z 1 ,z 2 ), 



(45a) 
(45b) 



5(3 2 (z 1 ,z 2 ) 

Eq. (45b) is the Dyson equation with suitably defined mass operator. 

In what follows it is supposed that the system under consideration is spatially ho- 
mogeneous, so that the correlators Pi(z) = ((f)) and f3 2 (z,z') = G(z,z') do not 
depend on the site index. Moreover, since Eq. (45a) has a trivial solution ((f)) = at 
Ai — and ergodicity breaking transition is related to the anomaly in G when FDT 
is violated, we eliminate ((f)) from the consideration by putting ((f)) =0. Eqs. (41- 
44) then provide the expression for the dynamical action: 



2$ = Trln[G] + dz 1 dz 2 G 4 (z 1 , z 2 ) + 



+ T" 



1 J &Zi&Z 2 



w 
2T 



G 2 (2 1? z 2 ) - K(z ± , z 2 )G(zi, z 2 ) 



(46) 
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where w — J2i The corresponding Dyson equation follows from the stationar- 

■• rt'I> 
lty condition — — = : 

dG(z, z') 

\ 111 

G-\z, z') = -K(z, z') - —G(z, z') - £ (z, z') . (47) 



The mass operator £ (z, z') = A 2 / (6T 2 )G 3 (z, z') has the following components 

E(f , = S(C) = ^r 2 C\t, f ) , E±(t, = £'(C) G ± , (48) 

, as 

where £' = — . Clearly, £ satisfies the FDT: d t ,£ = -d t £ = £_. 

By using the identity (22), it is not difficult to derive the components of Eq. (47). So 
the Dyson equation can be represented as a system of nonlinear integro-differential 
equations: 

Gl 1 = w 8(t- (m =F ft) - - £ ± , (49) 



GZ 1 ■ C ■ G- 1 = fa - + ^ + + £ . (50) 
Recall that the dot "■" denotes the operator product. 

Finally, we arrive at the explicit form of the dynamical Dyson equations for the 
autocorrelator C(t, t') and the response (Green) function G(t, t')\ 



±(m + d t )G(t, f ) = 5(t - + ^ \l dr' G(t, t')G(t' , + 

+ /*dr'E_(*,r')G(r',0 (51) 



i(m + ft)C(f , = |;G(f , t) + ^ dr' G(f , f) + 

f dr'C(t,r')G(t',r')+ f dr' G(t,r')C(r' ,t') 

J—oo J—OO 

+ [ l dr'j:(t,r')G(t',T r )+ f dT'Z_(t,T')C(T',t') 

J—oo J—oo 



W 
t' 



+ 



(52) 
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4 Ergodicity breaking transition 



In this section we concentrate on asymptotic analysis of the Dyson equations (51,52) 
Our primary purpose is to describe the behavior of the system in terms of static 
susceptibilities and memory parameters characterizing asymptotics of the autocor- 
relator C(t, t') in the limit of large time separation, r = t — t' — > oo. 



The first step is to transform Eq. (52) in such a way that FDT is not broken explic- 
itly. To this end we assume that interaction is adiabatically switched (C(— oo)=0) 
and apply FDT on the interval of time from — oo to t'. The result reads 



i (m + d t )C(t, if) = q { ^ + ^C(t, if) + E(f , if) | + 



+ Y2 1 dr ' G <t> W> + 1 dr ' s - (*» W> > (53) 



where go is the value of equal time autocorrelator C(t', t'). Thus the starting point 
for subsequent analysis is the Dyson equations (51) and (53). 



4.1 High temperature region 



Since in this region the time translational invariance is unbroken, Eqs. (51,53) can 
be easily analyzed by making use of Laplace transformation. The resulting system 
for the Laplace transforms C(p) and G(p) is given by 



-(m+p)G(p) = 1 + T2 G\p) + £_(p)G(p) , (54a) 

1 - ( h 2 w - « \ 

-[(m + p)C(p) - g ] = go f ^ + j^Cip) + £(p)J + 

+ C'( P )(^G( P ) + E_(p)) . (54b) 



By using the relations: \ — G(0), q = \im p ^ 00 pC(p), q h = lim p ^ pC(p) and 
£(g ft ) = \im p ^ pT l (p), from Eq. (54a) and Eq. (54b) multiplied by p the following 
system for susceptibilities and the field induced parameter can be deduced 
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m 


= 1 + 


W 


TjiX 


J^2 


m 




W 


TpXa 


= 1 + 


2^2 






W 






J>2 



Xa + Xa {S(Xa + Qh) ~ ^Qh)} ~ 



(0)} , 



(55a) 



(55b) 



[h 2 w } 

qh = x{Xa + qh) <7j^ + Tj^qh + ^iqh)} , (55c) 

where X a = <?o - q h , A X = X ~ Xa and A£_(0) = E_(0) + 5%) - S( Xa + 
qh). It should be stressed that, except for the identity lim p _ >0O p(pC'(p) — go) = 
— lirrip^oo pG{p) = — T that yield the first term on the right hand side of Eq. (55b), 
FDT have not been used in the derivation of the above equations. 

In the FDT regime Ax = and A£_ (0) = 0. So the equations for the susceptibility 
X and the field induced parameter q h in the high temperature region are 



-X = 1 + t^X 2 + X + qh) ~ Z(q h )} , 

qh=xix + qh) 1^2 + + • 



(56a) 
(56b) 



Notice that, in the case of h 2 = 0, the temperature dependence of x is defined by 
Eq. (56a) with q h = 0. This is why the parameter q h is referred to as a field induced 
parameter (there are no memory effects in the absence of random field at T > T c ). 



4.2 Critical point 



Ergodicity breaking transition is determined by the point where the FDT compliant 
solution with Cfdt(^) and GWr(£) becomes dynamically unstable. In order to 
study stability let us first differentiate Eq. (53) with respect to time and then insert 
the perturbed solution d t C = c^Cfdt + SD and G = GVdt + 5G into the resulting 
system. As far as stability analysis is concerned only linear part of the system is 
relevant. After rather straightforward calculations, it can be derived in the form: 



fdt + x' 1 - X + E'fah)) 



w 



5G = 
5D- 



-<„. [ — + ^(g h ))6G: 



(57a) 



(57b) 



where "... " stands for the terms that are nonlinear in 5D and SG. 
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As a result, the marginal stability condition reads 

l=x(x + 5h){^ + S'(?h)} • (58) 

It should be emphasized that fluctuations do not need to obey FDT, 5D + 5G ^ 0. 
As a consequence, at q h ^ Eq. (58) is stronger than the condition 

l = X 2 {^ + S'(g,)} . (59) 

obtained under the assumption that fluctuations do not violate FDT (SD + 5G = 0). 

Eqs. (56,58) yield the temperature of ergodicity breaking transition, T c , and the 
value of q h at the critical point, q c . The latter can be easily obtained from Eqs. (56b, 58): 

q c = (3/i 2 /A 2 )V3. 

4.3 Behavior at the limit T — > T c — 

Despite the detailed analysis of low temperature region is beyond the scope of 
this paper it is instructive to see how equation that define the limiting value of 
susceptibility when approaching T c from below is related to the marginal stability 
condition. 

Let us first consider Eqs. (55) and suppose that near T c FDT is modified as follows: 

~ C(t) = G(t) - AG(t) , AG(t) = -^(C(t)) , (60) 

where ^ is a nondecreasing function which derivative ^' vanishes outside the in- 
terval [91,92] with 92 < 9o- Obviously, these assumptions ensure the validity of 
Eq. (55) and lead to the following results: 

A X = tt(g 2 ) - , (61) 

AE_(0)=/ dgi;'(g)*'(g) = E , (g m )Ax, (62) 

Jqi 

where q m G (91, 92) is the middle point. 

Below T c Ax ^ 0, so subtracting Eq. (55b) from Eq. (55a) and dividing the result 
by Ax gives 

l = x(x + 9.){^ + S'(9 m )} , (63) 
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where the terms proportional to Ax are neglected. From Eqs. (55c, 63) q m — > qh 
at T — > T c — and we recover the marginal stability condition as equation for the 
limiting value of susceptibilities. 

Clearly, aging cannot be taken into account provided the time translational invari- 
ance is unbroken, so the above consideration is closely related to the regime known 
as true ergodicity breaking [22,23]. In this case the system equilibrates in a sepa- 
rate ergodic component. Alternatively, according to the concept of weak ergodicity 
breaking [24], the system does not equilibrate and asymptotically (t, t' 3> 1) FDT 
is violated at large separation times, r ~ t', whereas the system reveals quasi- 
equilibrium behavior at sufficiently small r. It implies that the generalized form of 
FDT breaking contributions to the correlation functions is [17] 

C(t, f) = C (t - f) + C a (t'/t), G(t, f) = G {t - t>) + t^Gait'/t), (64) 

and the modification of FDT is given by the following relations 

dtfCoit - = G (t - t'), xd z C a (z) = G a (z), (65) 

where z = t'/t and x parameterizes the violation of FDT. The time-dependent 
susceptibility x(M') (see Eq. (26)) now depends on both the waiting time t' and 
the separation time r. So the value of static susceptibility is different depending on 
the order in which limits t' — > oo and r — > oo are taken 

lim Jim x(t' + t, t') = Xa = qo-q, 

lim X (t' + r,t') = X = Xa + xAq, (66) 

T—fOO 

where q = C(t, t) = C o (0) + C (l), Aq = q - q h = C a {\) - C a (0) and 

q = lim lim C(t' + r, t') = C (oo) + C a {\) , 

r— i-oo f— ►oo 

q h = lim C(t' + r, t') = Co(oo) + C Q (0) . (67) 

In the sense of Cubo [25], Xa and x can be referred to as adiabatic (thermodynamic) 
and isothermal susceptibilities, respectively; q is the dynamical Edwards-Anderson 
parameter. Inserting Eq. (64) into Eq. (51) integrated over the second argument of 
G gives the system 



m 


1 + 


w 






m 


1 + 


w 


~7jiX 


J>2 



+ Xa{S( X a + g)-S(g)} , (68a) 
+■ X {nXa + q)- + {x- l)AqE'(q h )} . (68b) 

In exactly the same manner as the marginal stability condition (58) was recovered 
from Eqs. (55,61,62), Eqs. (68) yield the condition (59) in the limit Aq — > at 
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T — > T c — 0. Interestingly, it follows that, despite the susceptibility x an d non- 
ergodicity parameter Aq are continuous at the critical point, q can be discontinuous. 

So, in both cases A% plays the role of order parameter, but, in general, equations 
for the limiting values of susceptibilities are different depending on whether the 
aging is taken into consideration. Assuming the time homogeneity implies that the 
system can be described by Eqs.(55) below T c , but it can be shown that in this case 
we would encounter discontinuity of the susceptibility at the critical point. Thus 
we arrive at the conclusion that aging plays an important part in the problem under 
consideration. 



5 Numerical results and discussion 



In the above section the ergodicity breaking transition is treated on the basis of the 
dynamical approach. The critical temperature is defined as the point of marginal 
stability for FDT complaint solution. At h 2 = Eqs. (56a,58) are easy to solve, so 
the critical temperature T c is given by 



T c = 4 (l-Jl-w c /w), w c = 3 

fj,w c V v / 



2A 
3/x 



(69) 



From Eq. (69) it is clear that, in order for the transition to occur, the disorder inten- 
sity w must exceed its critical w c . Dependence of w c on the intensity of random field 
is depicted in Fig. 1. It is seen that there is no threshold for w at sufficiently large 
intensities of random field, whereas w c is an increasing function of h 2 in the range 
of small intensities. Note that, in order to simplify analysis, the relevant quantities 
are rescaled as follows: 



M, 



w 



2\-2 



6/TA 



w. 



h 2 -> 6/x 3 A~V T -> VQ^X^T , 



so the system for T c and Xc at the critical point takes the form: 



( X c - 1) T 2 = w X l + Xc (Xc + Qcf - q, 



T 2 c =Xc{Xc + q c ) w + 3q t 



(70a) 
(70b) 



where q c = h 2 ^ 3 /2. Dependencies of T c and the temperature of freezing transition, 
where the derivative of x with respect to temperature diverges, 7} on the disorder 
intensity w at h 2 = are shown in Fig. 2. As is seen, the difference T c — Tf goes to 
zero when w increases. The latter is not very surprising, for T c = Tf at A = and 
the rescaled intensity w increases indefinitely when A — > 0. As is shown in Fig. 3, 
dependence of T c on q c , that is proportional to h 2 / 3 , reveals more complicated non- 
monotonous behavior at w just above its critical value. The critical temperature 
becomes an increasing function of the random field intensity at sufficiently large w. 
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In solving the equations (56,58) one has to handle nonuniqueness of the solutions. 
For example, changing the sign of the radical in the expression for the critical tem- 
perature at h 2 = (69) gives another solution that corresponds to the nonphysical 
branch and has wrong behavior, T c — > oo, in the limit A — > 0. Note that choos- 
ing this branch gives another physically absurd result that q h grows as temperature 
increases at T > T c . 

Let us summarize the results of the paper. It is shown that the dynamical action of 
thermodynamic system with quenched disorder can be calculated as a second Leg- 
endre transformation of the effective free energy functional. Despite the technique 
was employed to study the disordered system in the mean field approximation, it 
can be equally applied for consideration of fluctuational effects when going beyond 
the scope of the mean field theory. From the other hand, owing to the reflection 
symmetry of the system, — > — 0, we have eliminated (0) from the considerations 
concerning the ergodicity breaking transitions in the case of symmetric distribution 
of quenched variables. Clearly, when the distribution is nonsymmetric or A 3 ^ 
(the presence of cubic anharmonicity), (0) ^ and the dynamical action contains 
the terms dependent on the averaged order parameter. Eqs. (45) then give coupled 
equations for the order parameter and correlation functions. 

The method is applied for the study of 4 model of thermodynamic system with 
quenched couplings and external field written in the site representation. Asymp- 
totics of correlator is found to be affected by the random static field. It is character- 
ized by the field induced parameter q h which is a decreasing function of temperature 
at T > T c . Discussion at the end of Sec. 4 led us to the conclusion that the differ- 
ence between the memory parameter q and qh plays the role of an order parameter, 
so that adiabatic and isothermal susceptibilities differ if Aq = q — q h ^ 0. Numeri- 
cal analysis reveals that using the system (55) below T c would predict discontinuity 
of the susceptibility at the critical point. Much more reasonable results can be ob- 
tained in the case of weak ergodicity breaking. So aging plays an important part 
in dynamics of the system below critical point. Recently reported results [27] for 
heteropolymers support the conclusion. 
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Appendix A 



In the bulk of the paper the solutions under investigation are casual. In general, this 
is not the case. An instanton motion, where both initial and final boundary condi- 
tions need to be fixed, provides an important example of this kind. In this Appendix 
we show that the mapping, recently introduced in [26], between the uphill motion 
related to an instanton and the corresponding downhill motion going back in time 
can be constructed by making use of SUSY formalism. 

Let us define the transformation of superfields as follows: 

<j>(z) - 0(z) = T_{z)^{z) , fa) - <f>(z) = T + (z)4(z) , (A.l) 

where 

T±(z) = exp(±ee fy) . (A.2) 

Inserting Eq. (A.l) into the action (9,10) integrated from t = U to t = tf gives the 
expression in terms of superfields <p 

S = S + (V f -Vi)/T, L = D_4>D_4> + V(4>), (A.3a) 
D_=T_{z)DT + {z) = - + 9- = D' M , (A.3b) 

D_ =T_(z)DT + (z) = t = D' M , (A.3c) 

where Vi = V(r](ti)) and Vf = V(r](tf)). From Eq. (A.l) it is clear that 

G( Zl , z 2 ) = T + { Zl )T + {z 2 )G{ Zl , z 2 ) , G( Zl , z 2 ) = . (A.4) 

The latter follows the relations: 



o2 c\ C) 

(V(*i)^2)> = Q^p(t x , t 2 ) + QfGi^h) + —G(t 2 , t ± ) , (A.5a) 

G(h, t 2 ) = G(h,t 2 ) + ^C(h, t 2 ) , (A.5b) 

ot 2 

C(t 1 ,t 2 )=C(t 1 ,t 2 ). (A.5c) 

It remains to notice that FDT is invariant under the action of the transformation 
(A.l) followed by the inversion of time t — > —t, so the Lagrangian L describes 
normal downhill motion formally inverted in time. The resulting relations (A. 5) are 
identical with those derived in [26] and enable calculating of the Green function for 
the instanton process provided that the corresponding causal solutions are known. 
More details on the subject will be published elsewhere. 
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FIGURE CAPTIONS 



Figure 1. Dependence of the critical value of quenched disorder intensity w c on 
the random field q c = [h 2 /2] 1 ^. (w c and h 2 are calculated in units of A 2 / (6/z 2 ) 
and A 2 / (6yU 3 ), respectively.) 

Figure 2. Dependencies of the critical temperatures T c (ergodicity breaking tran- 
sition) and Tf (freezing transition) on the quenched disorder intensity w in the 
absence of random field, h 2 = 0. (w and the temperatures are calculated in units 
of A 2 / (QfJ 2 ) and fj, 2 / {VQ\), respectively.) 

Figure 3. Temperature of ergodicity breaking transition T c versus the intensity of 
random field q c = [h 2 /2] 1 ^ at various values of w. 
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FIGURE 1. 
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FIGURE 2. 
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